Effect of Randomness on Quantum Data Buses of Heisenberg Spin Chains 
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A strongly coupled spin chain can mediate long-distance effective couplings or entanglement be- 
tween remote qubits, and can be used as a quantum data bus. We study how the fidelity of a 
spin-1/2 Heisenberg chain as a spin bus is affected by static random exchange couplings and mag- 
netic fields. We find that, while non-uniform exchange couplings preserve the isotropy of the qubit 
effective couplings, they cause the energy levels, the eigenstates, and the magnitude of the couplings 
to vary locally. On the other hand, random local magnetic fields lead to an avoided level crossing 
for the bus ground state manifold, and cause the effective qubit couplings to be anisotropic. Inter- 
estingly, the total magnetic moment of the ground state of an odd-size bus may not be parallel to 
the average magnetic field. Its alignment depends on both the direction of the average field and the 
field distribution, in contrast with the ground state of a single spin which always aligns with the 
applied magnetic field to minimize the Zeeman energy. Lastly, we calculate sensitivities of the spin 
bus to such local variations, which are potentially useful for evaluating decoherence when dynamical 
fluctuations in the exchange coupling or magnetic field are considered. 

PACS numbers: 03.67.-a, 75.10.Jm, 75.10.Pq, 75.75.-c 



-i— > 

a 

CM 
> 
00 
00 
(N 

o 

(N 



I. INTRODUCTION 

A qubit is the elementary unit of quantum informa- 
tion, and can be realized with a variety of two-level sys- 
tems, such as confined electron spins in a semiconductor 
nanostructure. For electron spin based qubits, universal 
quantum gates can be realized using Zeeman coupling 
and spin-spin exchange interaction.— The direct Heisen- 
berg exchange coupling between two electron spins is de- 
termined by the overlap of electron orbitals, and is thus 
a short-range nearest neighbor interaction. In order to 
implement quantum algorithms efficiently, quantum gate 
operations on remote qubits, i.e., controllable long-range 
couplings, are needed. Various quantum data buses have 
been introduced to bridge this gap£~— In this context, 
we have proposed to use the ground states of a strongly 
coupled spin chain as a quantum data bus, or a spin 
bus£ We have shown that the parity of the spin bus can 
significantly alter the long-range effective couplings and 
entanglement between qubits that are coupled to the spin 
bus;£ and external fields can modify the form of the ef- 
fective interaction between the attached qubits»i More 
recently, we have also shown that high-fidelity quantum 
state transfer can be achieved via such a spin bus<2 

An ideal quantum information processor has identi- 
cal qubits, with precise control over couplings between 
qubits, and the qubits should be well isolated from their 
environment. However, in reality it is essentially impos- 
sible to create identical qubits based on artificial struc- 
tures such as quantum dots and Joscphson junctions, and 
in a solid state environment there are normally several 
sources of qubit variance. For example, the size of a 
quantum dot and the electron orbitals are largely deter- 
mined by the gate structure and the applied gate volt- 
ages. They can also be strongly influenced by factors 
such as the band structure of the host semiconductor 



and the random potential landscape due to modulation 
doping. Furthermore, the Coulomb exchange coupling 
between spin qubits is determined by the exponentially 
small overlap of the electron orbitals, and controlled by 
the gate voltages. Small variations in gate voltages could 
thus cause large changes in the exchange coupling. Such 
deviations from the ideal value could lead to imperfect 
gate operations, and possible gate errors^ There are gen- 
erally also very slow charge traps in a semiconductor het- 
erostructure, where a trap can switch between two dif- 
ferent charge distributions at a time scale much longer 
than the qubit operation time scales. While such a trap 
would probably be static during a quantum operation, 
it could modify the exchange coupling to a value that is 
different from the calibrated value. Similarly, via hyper- 
fine interaction, environmental nuclear spins produce a 
local random magnetic field for a quantum dot confined 
electron spin qubit This field can be considered quasi- 
static in the context of a spin bus because its dynamics 
is much slower than the bus mediated gates. In short, 
in building a practical quantum information processor, 
deviations from calibrated values for various control pa- 
rameters are inevitable. It is thus necessary to know the 
engineering tolerance in the variation of parameters such 
as the spin-spin coupling and external magnetic fields. 

In this paper we study how the capabilities of a spin 
bus are affected by static random variations in the ex- 
change couplings between the bus node spins and the 
external magnetic fields experienced by the bus nodes. 
Specifically, we study how the bus spectrum, bus-qubit 
coupling, and bus-mediated qubit-qubit coupling are af- 
fected by these random but static variations of the system 
parameters. The paper is organized as follows. In Sec. |H] 
we discuss how the strongly coupled Heisenberg chain 
can act as a spin bus when external qubits are weakly 
attached to it. We derive the effective Hamiltonians of 
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the qubit-bus system up to second order. In Sec. IIII1 we 
show how the fluctuations in exchange couplings and ex- 
ternal magnetic fields within the spin bus could affect the 
fidelity of the bus. Finally, the summary and discussion 
are given in Sec. IIVI In the Appendices we discuss how 
to obtain the effective Hamiltonians using a projection 
method, and present more detailed results on the bus 
spectrum. 



II. SPIN CHAINS AS QUANTUM DATA BUSES: 
THE QUBIT-BUS EFFECTIVE HAMILTONIANS 

In this section, we discuss how a strongly coupled uni- 
form antiferromagnetic Heisenberg spin chain can be used 
as a quantum data bus, or spin bus, which coherently 
connects remote qubits^— In Appendix [Al using a many- 
body perturbation method based on the projection oper- 
ator, we derive the effective Hamiltonians for the qubit- 
qubit and qubit-bus couplings to first and second order. 
We also calculate numerically relevant energy gaps, local 
magnetic moments, and the effective couplings for finite 
spin chains. 
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FIG. 1: (Color online) Left panel: Two external qubits A and 
B (red circles) are weakly coupled to an (a) odd-size or (b) 
even-size spin bus (blue circles). The right-hand panels show 
the corresponding effective Hamiltonians in the low energy 
limit. An odd-size spin bus in its ground state acts as an 
effective spin- 1/2 particle denoted by C. It is coupled to the 
external qubits with strength j£) at first order in the per- 
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turbation theory, and induces an RKKY-like coupling J^J 
between the external qubits at second order. The even-size 
spin bus only mediates an RKKY-like coupling between the 
external qubits at second order. 

The system we consider has two spin qubits A and B 
weakly attached to an open spin chain C, as illustrated 
in Fig. [TJ The Hamiltonian of the total system^— is 



H = Hr 



H, 



QC 



Hr 



(2.1) 



In the ideal case with uniform exchange couplings and 
a uniform external magnetic field (along the z direction, 
B = B z), the Heisenberg Hamiltonian He of the spin- 
1/2 chain is written as 



JV-l 



H c = Jo Si ■ Sj 



N 



(2.2) 



where Sj is the spin operator of the i-th. node of the chain, 
N is the total number of spins in the chain, and Jo > 
indicates a uniform antiferromagnetic coupling between 
any two nearest neighbor spins at sites i and i + 1 . In this 
work, we will assume that a small or vanishing external 
magnetic field Bq is applied to the spin chain (The op- 
eration of a spin chain under a finite external magnetic 
field is discussed in Refs. and 123.). Hereafter, the spin 
chain C will be referred to as a spin bus because it acts 
as a quantum data bus. 

The Hamiltonian Hqc describes antiferromagnetic 
couplings of qubits A and B to the i-th spin and the 
j-th spin of the chain, respectively 



Hqc = Ja,% Si ■ Sa + Jbj Sj ■ S/ 



(2.3) 



We assume that the bus-qubit couplings J a j with a = 
A, B are small enough that the spin bus remains in its 
ground state manifold at all times. In this limit, the 
total Hamiltonian H can be split into the unperturbed 
Hamiltonian Hq = He + Hq and the perturbation Hi = 
Hqc- The perturbation condition is 



J Q ,,/A ~ N Jcm/tt 2 Jo < 1 



(2.4) 



where A ~ tt 2 J /N is the zero- field gap2£ above the 
ground state manifold (for an odd-size bus) or ground 
state (for an even-size bus). J Qj i/A is thus used as a per- 
turbation parameter. The qubit-bus coupling J a ^ can 
be turned on and off (gradually), unlike the static intra- 
bus coupling Jo. In general, external magnetic fields 
and B b may be applied to qubits A and B to implement 
single-qubit operations on them, so that the Hamiltonian 
Hq can be written as 



Hq = -gfJ-B^A ■ Sa - g^B^B • Si 



(2.5) 



There is no direct exchange coupling between qubits 
A and £?, as they are nominally well separated. As 
shown later, the spin bus C can mediate an effective 
coupling between them if they are both coupled to the 
bus. Since single-qubit operations are generally done 
separately from two- or multi-qubit operations, we set 
Hq = throughout this paper and will focus on qubit- 
bus and two-qubit couplings. Note that we will set ft = 1 
and Jo = 1 below for convenience. 

In the perturbative limit described by Eq. (|2.4|) , we can 
perform a canonical transformation of the full Hamilto- 
nian in Eq. (|2.1|) to obtain an effective Hamiltonian where 
the spin bus is in its ground state manifold. Details of 
the transformation are provided in Appendix A, as well 
as Refs. and [||. The actual form of the effective in- 
teraction depends on the parity of the bus. An odd-size 
spin bus has a doubly degenerate ground manifold, and 
acts as an effective spin-1/2 particle. At first order in the 
perturbation theory, the spin-1/2 bus couples directly to 
the external qubits. At second order, the spin bus medi- 
ates an effective RKKY-like coupling between the qubits. 
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The resulting effective Hamiltonian is given as follows^ — the perturbation is given by£ 



H 



(2) 
off 



jjg S B ■ S c + Jtf S A ■ S B 



r(2) 



(2.6) 



where Sc is a spin operator representing the ground dou- 
blet states of the spin bus. The effective coupling between 
qubit a and the spin-bus C is given to first order in the 
perturbation parameter by 



J (1) - J m 



(2.7a) 



= (0;i|a«|0;i) = -(0;-iK|0;-i) 



(0;ski 



10; 



(2.7b) 



This is a product of the bare coupling J a ^ between the 
i-th spin of the spin bus and the external qubit a and the 
expectation value m* of <Ji Z at site i in the ground state of 
the spin bus. Notice that although rrii is dimensionless, it 
can be considered as the local magnetic moment at site i 
when multiplied by g/is /2. The RKKY-likc second-order 
coupling j( 2 ' is given byii 



J 



(2) 



JaaJb, 



\n)(n\a jfl \0) 



Eq — E„ 



(2.8) 



Here E n and \n) are the eigenenergies and eigenstates of 
He of an isolated spin bus, and cr^ with \i = x,y, z stand 
for Pauli operators of the i-th spin of the spin chain. The 
prime symbol on the summation indicates the exclusion 
of the ground states. At zero external field, the ground 
state |0) in Eq. (f278|l can be either |0; \) or |0;-|), or 
any linear combination between them. The choice does 

{ 2) 

not change the value of J\j - At a finite magnetic field, 
however, if the ground state of the spin bus is degenerate, 
the two states are generally not spin- flipped image of each 
other. In this case, Eq. (|2.8|) has to be modified. 

The effective Hamiltonian (|2.6| shows that the odd- 
size bus at zero or low field acts as an effective spin-1/2 
particle that is coupled to the external qubits A and B, 
as illustrated in Fig. [TJ Although in general J^] » J- 2 j , 
the second order term plays an essential role in long-time 
evolutions, such as in quantum state transfer £ Thus our 
calculations in the rest of this paper are mostly concerned 
with these two coupling strengths. Furthermore, we focus 
on their normalized form )H t = / ' J a ,i and K^j = 

(2) 

Ji j Jo/ JaaJb, j, which depend only on the size N of the 
spin bus and the external magnetic field. 

For an even-size bus, the sub-Hilbcrt space of interest 
is spanned by the non-degenerate ground state of the bus 
and the four eigenstates of the two qubits (again we focus 
on the low-field limit). Within this space the bus does 
not have any dynamics as it is represented by a single 
ground state. As for the two qubits, there is no first 
order effective coupling between them here, in contrast 
to the case of an odd-size bus. The second-order qubit 
coupling term is obtained in the same way as for an odd- 
size bus. The effective Hamiltonian to second order in 



H 



(2) 



Jij- Sa ■ S_B 



(2.9) 



(2) 

where the RKKY-like coupling J> ■' has the same form as 
Eq. (|2.8p . In this case, the prime indicates that the non- 
degenerate ground state is excluded from the summation. 

Based on the effective Hamiltonians and the corre- 
sponding parameters, we can make some qualitative ob- 
servations on where the bus-qubit system might be sus- 
ceptible to randomness and fluctuations. In the case of 
an odd-size spin bus with attached qubits, the key fea- 
tures that determine the operation of the bus include the 
Zeeman splitting of the ground doublet of the spin bus, 
and the energy gap separating the ground doublet and 
the excited states, A12. The former depends on the mag- 
netic environment for the bus, while the latter depends 
on the interaction strength between the bus nodes. Both 

the qubit-bus couplings and the effective qubit-qubit 

(2) 

couplings ; • depend on the local exchange couplings 
and the local magnetic moments of the bus in its ground 
state manifold, which is a function of both magnetic en- 
vironment and the intra-bus exchange couplings. In the 

(2) 

case of an even-bus with attached qubits, J\ ? has simi- 
lar dependence on system environment as in the odd-size 
bus case, and is thus susceptible to variations in both the 
local magnetic fields and exchange couplings. 



III. EFFECTS OF RANDOMNESS 

In SecHH we have shown how a Heiscnbcrg spin chain 
with uniform exchange coupling Jo acts as a spin bus. 
Now we address the main question of the present paper, 
on how static randomness in exchange couplings and ex- 
ternal magnetic fields can affect the fitness of the spin 
chain as a quantum data bus. More specifically, we in- 
vestigate how such randomness influences the two energy 
gaps, Api and A12, and the effective qubit-bus and qubit- 



qubit couplings J}*] and j[ 2 ' 



In order to take into account the effects of randomness 
in exchange couplings and applied magnetic fields, the 
Hamiltonian of the chain (|2.2p is generalized to 



H c = 



N-l 

E 



Ji Sj • Sj + i - g^ B 



N 

E 



Bi 



(3.1) 



where Ji > is the antifcrromagnctic coupling between 
two neighboring spins at sites i and i + 1, and Bi is the 
local magnetic field at the ith site of the spin bus. Note 
that in spite of the random Jj and 73j, it can be easily 
shown that Hamiltonian ( 13. ip still commutes with the z 
component of the total spin S, i.e., [He, S z ] = 0. 

Hamiltonian Q3.ip may be considered as a finite quan- 
tum spin glass model.— There are several spin glass mod- 
els depending on the types of couplings (Ising or Hcisen- 
berg, and short range or long range) and the distribu- 
tions of Ji. For example, the Shcrrington-Kirkpatrick 
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models has the couplings between arbitrary pairs, sam- 
pled from the normal distribution with zero mean, while 
the Edwards-Anderson modet^ has only the nearest 
neighbor couplings. In the context of quantum infor- 
mation processing, we can reasonably assume that the 
exchange couplings and the applied magnetic fields are 
both near their target values, Jo and Bq. The random 
exchange couplings Ji and applied magnetic fields Bi are 
then 



(a) 



Ji 



J + 5J t , 
B Q + SB, 



(3.2a) 
(3.2b) 



In the numerical analysis described below, we choose SJi 
and SBi that are randomly sampled from normal dis- 
tributions with standard deviations a,; and erg, respec- 
tively. In experimental systems, we would expect such 
variations to be small, assuming reasonable calibration 
efforts. In the following studies, we analyze these two 
types of random variations separately, keeping one of the 
variables uniform. 



Effects of Random Exchange Couplings in 
Odd-size Buses 



In this subsection, we investigate how random varia- 
tions in the inter- node exchange couplings Ji, given by 
Eq. p. 2b ). affect the ability of an odd-size chain to func- 
tion as a spin bus. Such variations could result from 
calibration errors, slow but random hopping of charge 
traps near the spin bus nodes, and whatever other factors 
that are not accounted for during the calibration process. 
Here the external magnetic field Bo on the chain is set 
to be zero or small, so that the system remains in the 
isotropic regime. The external magnetic field, if any, is 
taken to be uniform, so that SBi = 0. 

The variations SJi are small compared to Jo, so that 
they may be treated as a perturbation: 



H C = H < § ) + V, 



(3.3a) 



where is the unperturbed Hamiltonian (|2.2p . and 

the perturbation V is given by 

N-l 



V = ^2 SJi St ■ s i+1 



(3.3b) 



Hereafter the superscript (°) is used to denote the case of 
uniform exchange coupling or uniform magnetic fields in 
the bus. 

Here we consider an ensemble of N$ odd-size chains. 
Each of the samples in this ensemble has the same size 
N but different J 2 ; sampled from the normal distribution 
with average Jo and standard deviation <tj. Fig. [2ja) 
shows the fluctuations in energy levels as a function of 
sample number, at a low magnetic field of gfisBo/ Jo = 
0.1. The two lowest energy levels Eq and E\ fluctuate in 
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FIG. 2: (Color online) (a) Four lowest energy levels, 
-Bo, Ei, E2 and £3, and (b) two energy gaps, Aoi = E\ — Eq 
and A12 = E2 — Ei are plotted as a function of sample num- 
ber Ns of odd-size chains with size N — 7 with random ex- 
change couplings Ji with the standard deviation aj/Jo = 0.1 
at gusBo/ J$ — 0.1. In (c) and (d), two energy gaps Aoi and 
A12 are plotted as a function of So for the samples N„ = 100. 



sync, so that the gap Aoi is free from the randomness of 
Ji [shown in Fig. [2(b)]. The energies E 2 and E 3 of |2) 
and 1 3) are also in sync, as shown in Fig. [5] (a). However, 
the gap A 12, which is a measure of the isolation of the 
ground doublet from the excited states, does fluctuate 
[as shown in Fig. |2(b)] , because E\ and E2 have different 
dependence on the exchange coupling. In other words, 
while the ground state splitting A i of this odd-size bus 
is robust against the randomness in exchange coupling, 
the ground-excited-state gap A12 is affected by the ran- 
domness. In Figs. [2] (c) and (d) the two gaps, Aoi and 
A12, arc plotted as a function of the uniform magnetic 
field Bq applied on the bus. At low fields, the ground 
state gap Aoi increases linearly and without broadening 
as the magnetic field increases, until g^isBo/ Jo ~ 0.35. 
It starts to be influenced by the exchange randomness 
above g/is-Bo/Jo ~ 0.35, which corresponds to the cross- 
ing between levels |2) and |3), as shown in Fig. [TJ] in 
Appendix IA II Beyond this crossing point, states |1) 
and |0) are not the time- reversal of each other anymore 
due to the level crossings with higher excited states. Re- 
call that for a spin chain to act as a spin bus, we need 
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the ground state doublet to be well separated from the 
excited states, or A12 3> A01. Panels (c) and (d) of 
Fig. [2] indicate that this condition is satisfied when the 
bus ground doublet is energetically separated from ex- 
cited states (gusBo/ Jq < 0.2), and acts as an effective 
spin- 1/2 system with a constant magnetic moment. 

Although [He, S z ) = dictates that the dimensionlcss 
total magnetic moment njj = ±1 of the ground state 
is still a good quantum number despite random exchange 
couplings, the local magnetic moment m, does fluctuate 



round 



as shown in Fig [3] Consequently, the first- 
order effective coupling between the qubit and the 
bus, given by Eq. (|2.7p . is affected by the randomness 
in the bus exchange coupling, and has to be calibrated 
individually. 
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FIG. 3: (Color online) Dimensionless local magnetic moments 
rrii = (0; ||o- 4z |0; |) are plotted (red dots) for all lattice sites 
on an odd-size spin bus with N = 9 for an ensemble of Ns = 
100 random samples, with ctj/Jq = 0.01 and Bo = 0. The 
height of the blue bars represents the local magnetic moment 
mf 1 for uniform exchange couplings (6 J = 0). The spread of 
the red dots at each bus node illustrates the variations in local 
magnetic moments due to the random exchange couplings, 
with larger variations apparent at the end sites. 



(2) 



is also affected 



The effective qubit-qubit coupling J> 
by the randomness of the intra-bus exchange couplings 
Ji. As indicated in Eq. (|2.8p . J-^ is determined by 
the full bus spectrum, including both the energy lev- 
els and the excited states. Figures [2] and [3] show that 
the random exchange couplings Ji in general affect the 
energy gaps from the ground state, E n — Eq, as well 

(2) 

as the bus eigenstates |n). Thus we expect that J!> ■' 
should be sensitive to the randomness in Ji. Figure 0] 
shows how the ensemble averages of the gap, (A i) e n, 
the local magnetic moment, (7715 ) on (which is the normal- 
ized first-order qubit-bus coupling J^\/ J A. 5), and the 
normalized second-order effective coupling, (Kij) cn = 

(Jij)en Jo/ (J A, iJb,5) depend on the fluctuations of the 
exchange coupling, represented by the standard deviation 
a j, over a 5-nodc bus. The blue filled circles in Fig. [4] 
show how the fluctuations in these quantities depend on 
the randomness in the exchange couplings. For example, 
when cfj/Jq = 0.1, 17(7715) ~ v{K\,s) ~ gj which indi- 
cates that the effective qubit-bus coupling and the effec- 
tive qubit-qubit coupling have similar sensitivities to the 
random variations in the intra-bus exchange couplings. 



Furthermore, both the local magnetic moment (thus 
the qubit-bus coupling) and the effective qubit coupling 
are linear functions of oj, with their slopes depending 
on the size of the bus. These slopes are indicators of 
sensitivity of and J^ to the exchange variations, 
and can be used in evaluating decoherence in such a spin 
bus architecture. For example, background charge fluc- 
tuations can affect exchange couplings between neigh- 
boring nodes of a spin bus. As a result, the effective 
qubit-qubit exchange coupling becomes a time-dependent 
random variable, which leads to two-qubit dephasing.— 
The relevant correlation function that determines the de- 
phasing is (j( 2 \t) j( 2 \0)), and is given approximately by 

[<j(JW)/<j,j] 2 (Ji(t)Ji(0))&. The latter correlation func- 
tion, (Ji(t) Jj(0)), represents fluctuations in the individ- 
ual inter-node exchange couplings along the bus, whose 
dynamics is determined by the environmental charge 
noise. 
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FIG. 4: (Color online) For odd-size chain with N = 5, the 
ensemble averages of (a) the gap (Aoi} cn in units of Jo, (b) 
the dimensionless local magnetic moment (ms) eD at the end 
of the chain, and (c) the normalized second-order effective 

coupling {Kij} cn = (J 1 <2 5) C nJo/(J'A,i^s,5) are plotted (open 
red circle) as a function of the standard deviation aj of SJi . 
The blue filled circles in each panel indicate the standard 
deviations of each data point, which is obtained by averaging 
over 2000 random configurations. Here Bo = 0. 

In summary, even when intra-bus exchange couplings 
of an odd-size spin chain have random but static varia- 
tions, the chain can still act as a spin bus, with a ground 
state doublet that is well separated from the excited 
states, and acts as an effective spin-1/2 system with a 
constant magnetic moment. However, the effective qubit- 
bus couplings and the mediated qubit-qubit couplings are 
affected by the randomness in exchange, with their flue- 
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tuations linearly proportional to the randomness in ex- 
change. Calibration would thus be needed for accurate 
qubit operations. The results here also have implications 
for spin-bus related decoherence. In essence, the strong 
exchange couplings allow a spin bus to process quantum 
information across a large distance, but also make the 
qubit-bus system susceptible to charge noise via both j\ ^ 
and jf}. 

B. Effects of Random Magnetic Fields in Odd-size 
Buses 

Spin qubits are generally susceptible to magnetic noise, 
and the spin bus is no exception. Here we examine 
how random but static external magnetic fields affect the 
properties of an odd-size spin bus. Our results should 
also be a useful indicator of the sensitivity of a spin bus 
to temporal magnetic noise, as we will discuss later in the 
section. For this calculation we assume that the exchange 
couplings Ji are uniform, and focus on the magnetic ran- 
domness. 

The local magnetic fields, Eq. (|3.2b[) . are 

B t =B + 5Bi , (3.4) 

where the random field SBi is sampled from a normal dis- 
tribution with standard deviation <tb- As in the previous 
subsection, we consider an ensemble of Ns spin chains, 
so that the ensemble average of local magnetic field Bi 
is (Bi) cn = B and (5Bi) cn = in the limit of large 
Ns- Such a random distribution of local magnetic field 
could originate from quasi-static nuclear hyperfine fields, 
or local paramagnetic centers in a semiconductor. 




9V B B o /J o 

FIG. 5: (Color online) Ensemble average of the gap (Aoi} en 
in units of Jo between the two lowest states of a single spin as 
a function of the external magnetic field Bo- The height and 
the width of the rectangle, touching the bottom of the curve, 
indicates g/J,B&B/Jo and (Aoi} on /</o at Bo = 0, respectively. 

As a benchmark, we first recall how an odd-size spin 
chain behaves in a constant uniform magnetic field. 
When Bq = 0, the odd-size chain has two doubly de- 
generate ground states with the total magnetic moment 
S z = ±1/2. A finite B splits these two states like a sin- 
gle spin- 1/2 particle. There are two important features 
that determine the behavior of the ground doublet in a 



magnetic field: the energy splitting Aoi, and the ground 
state spin orientation. In a uniform field the latter de- 
pends only on the <?-factor of the material, while in a 
random field it also depends on the local field configura- 
tion. 

We first examine how the random external magnetic 
field affects Aoi, the splitting of the ground state doublet 
of the odd-size bus. This dependence is a good indicator 
of whether the ground state doublet of an odd-size spin 
bus is a robust effective spin- 1/2. Note that the commu- 
tation relation [He, S z ] = holds even in a random mag- 
netic field, so the two lowest states have S z = ±1/2 (re- 
call that the inter-node coupling is anti- ferromagnetic) . 
Figure [5] shows that the randomness in the magnetic field 
induces a finite average gap at zero field. This gap opens 
at i?o = even though (SBi) cn = 0. This non- vanishing 
average gap between the two lowest states is a conse- 
quence of the statistical behavior of the random field, and 
can be understood using a model of Ns single spins in 
a Gaussian ensemble of random external magnetic fields 
8B with Bq = 0. The Hamiltonian of a single spin is 

H = - gf i B SB ^ . (3.5) 

The energy splitting of each spin is given by Aoi = 
gfx B \6B\, which is always positive. It is thus not a sur- 
prise that the ensemble average of the gaps, (A 01 ) cn = 
gHB(\SB\) on , is nonzero, even though (SB) cn = — the 
ground state changes according to the field configuration. 

To make this argument more rigorous, recall that for 
a normal distribution, the odd central absolute moments 
of a random variable X with a mean of /Lt are given by 

a p (l - «)!! 

E(\X-n\*)= I W , (3.6) 

V Z7T 

where !! denotes the double factorial. Applying Eq. (|3.6p 
to Aoi, we get 

(Aoi)cn = g/J.B-^== , (3.7) 

V 27T 

where as is the standard deviation of SB. While in our 
case the magnetic moment of the ground state is dis- 
tributed throughout the whole spin chain, its splitting 
is still mainly due to the Zeeman splitting of a single 
Bohr magneton, so that the single-spin argument pro- 
vided here is still applicable. When Bo / 0, the av- 
erage gap can be qualitatively expressed as (Aoi)en °t 
g\i B \J{(Bq + SB) 2 ) ~ gfi B yjB'l + a%, which is a hyper- 
bola that saturates at Bq = to g[iB<^B and approaches 
g/is£>o at large Bq. The difference in the proportionality 
constant originates from the difference between (|/|) and 

Now we address the spin orientation of the odd-size 
spin bus in the ground state under random magnetic 
fields. The magnetic moment of a single spin (with g = 2) 
in the ground state is anti-parallel to the external field. 
As in the single spin case, the total magnetic moment 
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FIG. 6: (Color online) (a) Total magnetic moment m of 
the ground state |0), (b) local magnetic moment mg for the 
ground state |0) (orange solid circle) and for the first ex- 
cited state |1) (green open circle), and (c) the fidelity F of 
the ground state in the presence of random magnetic fields 
with respect to that in the uniform field as a function of 
the sample number. Here the size of spin buses is N = 9, 
g^,BBo/Jo = 0.1, and g/iBO-B / Jo = 0.065. 
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FIG. 7: (Color online) The fraction of flipped ground states 
for N s = 10000 random realizations of odd-size chains with 
N = 5, as a function of Bo and as- 



of the odd-size spin bus in the ground state is always 
anti-parallel to a small uniform external magnetic field. 
As shown in Fig. [31 the local magnetic moments align 
antiferromagnetically (alternation between anti-parallel 
and parallel alignments). Even in a random external 
magnetic field in the z-direction, the two lowest states 
of an odd-size spin bus are characterized by S z = 1/2 or 



S z = —1/2, so one might guess the Zeeman energy would 
determine the spin orientation. However, this is not the 
case. 

Figure [H] shows the total magnetic moment of the 
ground state |0) in an ensemble of random magnetic 
fields for the case when N ~ 9 and <jb/Bq = 0.65. As 
illustrated in panels (a) and (c), seven samples, #10, 
#19, #36, #44, #57, #69, and #97, among 100 sam- 
ples have negative total magnetic moment m = — 1, so 
that their ground state is \n = 0; S z = — i) instead of 
\n = 0; S z = |). Panel (b) of Fig. [5] shows a strongly ran- 
dom local magnetic moment on the 9-node spin bus, with 
the seven flipped states having reversed local magnetic 
moments. The probability for either \n = 0; S z = — |) or 
| n = 0; S z = |) to be the ground state depends on the 
ratio of B and erg, as illustrated in Fig. [71 As expected, 
when Bq = 0, the probability is 50%. A strong uniform 
field Bq (compared to as) suppresses the flipping frac- 
tion and stabilizes one of the states as the predominant 
ground state. 

A close inspection of the data presented in Fig. [6] re- 
veals that Zeeman splitting of the local nodes of a spin 
bus does not tell the whole story of Aio, so that the 
single-spin model has its limitations. By definition, the 
energy gap is 
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(3.8) 



Here Hj and Hz represent the exchange and Zeeman 
components of H c in Eq. (|3~Tj) . E% = (k\H z \k) = 
-/i^^- Bi-m\ and Ej = (k\Hj\k) denote the Zeeman 
and exchange contributions to the energy of state |fc), 
with k = 0,1, = {k\<Ji Z \k) is the local magnetic mo- 
ment at bus node i, and we have taken g — 2 for simplic- 
ity. When the external magnetic field is uniform, Bi = 
Bq, the exchange contribution to the energy of |0) and |1) 
in Eq. (|3.8|) are the same, so that Ej — Ej = 0, and m- 
are equal in magnitude and opposite in sign to m\ . Thus 
the gap between the two states is given completely by the 
Zeeman splitting: E\ — = 2li>BBJ2i m i = ^^bBtvP . 
With g = 2, the net spin of the ground state |0) is paral- 
lel to the uniform magnetic field B like a single spin, so 
that |0; S z = 1/2) is the ground state when B is positive. 

Under a random magnetic field Bi, but with Bq 3> <tb, 
the exchange contribution to the gap in Eq. (|3.8|l is still 
negligible, and the ground state spin is parallel to Bq. 
However, if Bq is of the order of erg or smaller, the 
ground state total spin orientation may be anti-parallel 
rather than parallel to Bq . For these configurations with 
the "flipped" ground state, the reason for the flipping 
is varied. Consider the samples (from the ensemble of 
9-node buses presented in Fig. [6]) shown in Fig. [8] and 
Table [U Here sample #0 refers to the spin bus in a uni- 
form magnetic field. Samples #6, #10, #19 and #36 
are for the bus in different random magnetic field con- 
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figurations, and the later three samples have a flipped 
ground state, with m° = — 1. Figure [8] shows that in a 
random magnetic field, the local magnetic moments for 
the two lowest states are generally not equal in magni- 
tude, |m°| # l^l, so that the two states |0) and |1) arc 
no-longer spin-flipped image of each other, even though 
J2i m i = The most dramatic examples are those 

samples with a flipped ground state, where the ground 
and first excited states are far from the classical anti- 
ferromagnetic spin configurations. This also implies that 
in general the Zeeman energy E\ # —Eg. Samples #6 
and #10 demonstrate that even if all the magnetic fields 
Bt are positive, the ground state could still be either 
S z = 1/2 or S z = —1/2. While in most cases the energy 
gap between the two lowest states in Eq. (|3.8|) is due to 
the Zeeman term, for samples #10 and #19 the Zeeman 
energies E\ and E l z are almost equal, so that there is no 
Zeeman gap. The energy gaps for these two samples are 
determined by the exchange contribution in Eq. (|3.8|l . 
and are one or two orders of magnitude smaller than 
the usual Zeeman gap. For sample #36, the total en- 
ergy gap is dominated by the Zeeman contribution (with 
significant contribution also coming from the exchange 
interaction), although the configuration of the random 
magnetic field is such that the ground state is flipped. 
In short, when Bq < ob, the local and total magnetic 
moments of the ground state are sensitively dependent 
on the distribution of the random magnetic field. 



Sample No. 


Ei — Eo 


Ez — E° z 




\AEj/AE z \ 


#0 


0.2 


0.2 


0.0 


% 


#6 


0.27083 


0.25840 


0.01244 


4.8 % 


#10 


0.00403 


-0.00840 


0.01243 


148.0 % 


#19 


0.06222 


-0.01582 


0.07804 


493.3 % 


#36 


0.06989 


0.04306 


0.02683 


62.3 % 



TABLE I: (Color online) Energy gap Ei — Eo between the 
bus states |0) and |1), its Zeeman contribution E x z — E%, the 
exchange contribution Ej — E®, and the absolute ratio of 
these contributions for samples #0, #6, #10, #19, and #36 
in Fig. [6] Here the energy is measured in units of Jo. 

Our results so far indicate that random magnetic fields 
can seriously undermine the capabilities of a spin bus, by 
altering the local magnetic moments (and thus the effec- 
tive qubit-bus coupling jW) and the bus ground state. 
Fortunately, in general Ob is relatively small. For ex- 
ample, if the random field is due to hyperfine interac- 
tion in GaAs, ob ~ 2 mT for a 100 nm quantum dot, 
so that a Bq > 20 mT should be more than enough to 
overcome the effect of the random field. Furthermore, 
the above discussion is applicable to the regime where 
the magnetic energy scales are not much smaller than 
the exchange energy scales (e.g., for Fig. [6] and Table U 
HbOb/Jq = 0.0650). If Jo 3> B and ob, the structures 
of the ground and first excited states of the spin chain 
should be determined by the anti-ferromagnetic coupling 



and are less susceptible to the small magnetic field or its 
fluctuations. In this case the ground state spin orienta- 
tion would be mostly determined by the Zeeman contri- 
bution to the energies of the two states, and we would 
recover the simple single-spin physical picture. 
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FIG. 8: (Color online) Local magnetic moments of the ground 
and first excited states for several random field configurations 
in a 9-node spin bus. Here sample #0 corresponds to the 
case of a uniform applied field. The others correspond to the 
samples in Fig. [6] The total magnetic moment of the ground 
states of samples #0 and #6 is m° = 1, i.e., S z = 1/2. On the 
other hand, samples #10, #19, and #36 have the spin-flipped 
ground states, m° = —1, i.e., S z = —1/2. 

The fluctuations in the local magnetic moments of the 
spin bus due to the random external magnetic field lead 
directly to fluctuations in the effective qubit-bus cou- 
plings J^ [see Eq. (|2.7|) ]. though this fluctuation is sup- 
pressed if Bq 3> ob ! as illustrated by the few data points 
for small ob/Bq in Panel (b) of Fig. [§J In addition, un- 
der a random magnetic field, the effective coupling J^ 1 ], 
given by Eq. (|2.7p . becomes anisotropic. This is in addi- 
tion to the anisotropy induced by a finite -Bo-~ I n general, 
anisotropy occurs when 

1 



(0k«|l)# - 



(0|<7 iz |0> 



(lOiJl) 



and 



i=i 



|i>^i 



(3.9a) 



(3.9b) 



The anisotropy appears in both J^ and j( 2 )— , as illus- 
trated in Panels (b) and (c) of Fig. [Sj Based on the 
standard deviation data presented in panel (b), we also 
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observe that while the transverse component of the lo- 
cal magnetic moment is reasonably robust against the 
randomness in the magnetic field, the longitudinal com- 
ponent is not. On the other hand, panel (c) shows that 
both the longitudinal and transverse components of the 
effective qubit-qubit coupling J*- 2 ) have a linear depen- 
dence on the field randomness for small ub, changing to 
a different slope as <jb becomes larger than Bq. Both of 
these observations illustrate the fact that for an odd-size 
spin bus to function properly, magnetic field randomness 
in the system needs to be minimized. 



in the z-direction on an odd-size bus are as follows. First, 
the degeneracy in the ground states of an odd-size bus 
is lifted. Second, although the two lowest states have 
either S z = 1/2 or S z = —1/2, the spin orientation of 
the ground state is not solely determined by the Zee- 
man energy. Third, the random magnetic fields make 
the first-order and the second-order effective couplings 
anisotropic. 



C. Effects of Random Exchange Couplings in 
Even-size Buses 
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FIG. 9: (Color online) For an odd-size spin bus with N = 5, 
the ensemble averages of (a) the gap (Ai2) cn in units of Jo, 
(b) the local magnetic moment (rri5}en on one end of the 
chain, and (c) the normalized second-order effective coupling 
(if 1,5) en are plotted as a function of the standard deviation 
<jb of 8Bi. In (b), the labels "X" (symbol: red open cir- 
cle) and "Z" (symbol: orange open diamond) stand for the 
x- and z- components of the first-order effective coupling, re- 
spectively. In (c) the labels "XX" (symbol: red open circle) 
and "ZZ" (symbol: orange open diamond) refer to the 11- 
and 22-component of the second-order effective coupling, re- 
spectively. Here /xs-Bo/Jo = 0.05. In panels (b) and (c), the 
solid diamonds and solid circles represent standard deviations 
for the data labeled by the open diamonds and open circles, 
respectively. 

As explained in Sec. [Til the effective qubit-qubit cou- 
pling j( 2 ) given by Eq. (|2.8|) is well defined when the de- 
generate or nearly degenerate ground states of the spin 
bus are spin-flipped states of each other. This is the case 
for an odd-size spin bus near zero external magnetic field. 
At certain finite external fields, the ground states of the 
spin bus would be close to be degenerate again. However, 
in those regimes Eq. (j2.8|) is generally not applicable and 
should be modified^ 2 - 

In summary, the effects of the random magnetic field 



In this subsection, we investigate how an even-size spin 
bus is affected by random exchange couplings Jj. Recall 
that the effective Hamiltonian for an even-size bus cou- 
pled to two qubits in zero magnetic field takes the form of 
H cS = J S S a-Sb, where J$ is given by Eq. (23). The 
various terms in this equation are the bare qubit-bus cou- 
plings Ja,% and Jb.j, the bus energy gaps Ao„ = E n — Eq, 
and the transition matrix elements (0|<7j„|n). Here we 
examine the effects of random intra-bus exchange cou- 
plings on these different terms, with a particular focus 
on the gap A01 between the ground and the first excited 
state, which is an indicator of how well the nondegener- 
ate ground state is isolated from the excited states, and 
figures prominently in the expression of J*- 2 ' . 

An even-size bus can be thought of as an odd-size bus 
plus an extra spin- 1/2 node. Thus the lowest four states 
of an even-size bus is a singlet and a triplet. A uniform 
magnetic field would split the triplet, but would not af- 
fect the singlet ground state. With random exchange 
couplings, similar to the case of an odd-size bus, we still 
have [H, S z ] = 0, so that the triplet splitting is given 
by Zeeman splitting. The random exchange couplings 
do cause the energy levels and the eigenstates to vary in 
general. For example, Fig. 1101 plots the two energy gaps, 
A01 and A12, as a function of Bo for 100 samples of even- 
size chains with random exchange couplings. While the 
ground state remains nondegenerate, the gap A i is now 
distributed between 0.3 and 0.5 Jo at Bq = 0. Around 
Bq = 0, the gap A12 and the next gap A23 = -E3 — E% 
(not plotted) are robust against the random exchange 
couplings Jj, reflecting the fact that these two gaps cor- 
respond to the Zeeman splittings of the triplet bus states. 
This is similar to the A01 gap of an odd-size bus near zero 
field, which corresponds to the Zeeman splitting of the 
spin- 1/2 doublet ground state, as shown in Fig. [2] 

The most important effect of the randomness in the 
exchange couplings for bus operations is on the low ly- 
ing energy gaps such as A01. As shown in panel (b) of 
Fig. QTJ the ground state itself is quite robust against 
the random exchange coupling in terms of state fidelity, 
and has zero local magnetic moments as well as zero to- 
tal magnetic moment. This is in contrast with the effect 
of random magnetic fields on the even-size bus as shown 
in subsection IIII D[ where the local magnetic moments 
become non-zero. On the other hand, while the aver- 
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FIG. 10: (Color online) For N s — 100 even-size chains 
with random exchange couplings, two lowest energy gaps, (a) 
A12 = E2 — Ei and (b) A01 = E\ — Eq, are plotted as func- 
tions of Bo- The size of each chain is N — 8 and the standard 
deviation uj/Jq of SJi is 0.025. 

age value of the A01 gap is only weakly dependent on 
aj, the fluctuations in A01 depend linearly on aj, as 
shown in panel (a) of Fig. QT] Consequently, the fluctu- 
ations of the bus-mediated qubit-qubit coupling also has 
a linear dependence on the standard deviation aj of SJi, 
as shown in panel (c) of Fig. QTJ The slope <7jm I & J is 
quite large here, reflecting a sensitive dependence of the 
singlet-triplet splitting of the spin bus on the intra-bus 
exchange couplings. While calibration should be able to 
largely suppress the effects of any static randomness of 
the exchange coupling, the sensitivity to randomness in 
the local exchange, as indicated by the large aj(2)/crj, 
dictates that effects of environmental charge noise on the 
inter-node exchange coupling Ji need to be minimized. 



D. Effects of Random Magnetic Fields in Even-size 

Buses 

In this subsection, we study how an even-size bus is 
affected by an external magnetic field that has random 
local variations. In a uniform but small external mag- 
netic field, the ground state of the even-size bus has zero 
local magnetic moments as well as zero total magnetic 
moment (S z = 0). Recall that the commutation relation 
[He, S z ) = is still valid even when the bus is subject to 
random magnetic fields Bi in the z-direction. If the ran- 
dom magnetic fields are weak, the even-size bus remains 
in the ground state with S z — 0, i.e., zero total mag- 
netic moment. However, the local magnetic moments 
mi become non-zero in contrast with the case of a uni- 
form magnetic field. Furthermore, the bus excited states 
generally do have net magnetic moments, so that they 
respond to both local and global magnetic fields in terms 
of their energies and their state composition. 

In Fig. [T^] we plot the effect of the locally random mag- 
netic field on the gap Aqi, the ground state robustness 




FIG. 11: (Color online) For an even-size bus with N = 6, 

(a) the ensemble averages of the ground energy gap A01, 

(b) the fidelity F with respect to the uniform ground state, 
and (c) the normalized effective qubit-qubit coupling 7^1,6 = 
JyeJo/ J\,aJq,b are plotted as functions of the fluctuation aj 
of the random exchange coupling (all data are represented 
by red open circles). Here Bo = and the sample size is 
N s = 100. The blue solid circles represent the standard devi- 
ation of Aqi, F, and -R'i.e, respectively. 

(in terms of the fidelity F with respect to the uniform 
field ground state), and the bus- mediated qubit-qubit in- 

(2} 

tcraction J{ g . The average of gap A01 decreases with 

<jb because for any particular random field configuration, 

one of the polarized triplets has a lower energy compared 

to the zero field case and becomes the first excited state. 

The linear increase in the standard deviation of A01 is 

simply a reflection of the linear nature of the Zeeman 

splitting. In panel (b), the decrease in the fidelity shows 

that the local magnetic moments may fluctuate although 

their sum, i.e., the total magnetic moment is still zero. 

Interestingly, in panel (c), the normalized effective cou- 
(2) 

pling Ki.Q — J{ j [Jo/ Ja,iJb,6 and its standard deviation 
both have a super-linear dependence on a j when it is 
small, making the coupling a robust quantity against field 
randomness. The qualitative reason for this robustness is 
that the ground state is not coupled to the two polarized 
triplet states of the bus by the random magnetic field 
along z direction at the lowest order, while the gap be- 
tween the unpolarized triplet state and the ground state 
only depends on the field randomness quadratically. In 
this calculation we did not consider a random transverse 
field since our focus is on static disorder. For example, 
for magnetic disorder caused by random nuclear polariza- 
tions, the transverse polarizations would precess around 
the external field, so that their effect would tend to be 
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FIG. 12: (Color online) Effect of random magnetic field for an 
even-size bus with N — 6. (a) The average ground energy gap 
Aoi in units of Jo, (b) the average fidelity F of the ground 
states with respect to that without random Bj, and (c) the 
average of the normalized second-order effective coupling Ki,e 
are plotted (in open red circles or open orange diamonds) as 
a function of the standard deviation of the random magnetic 
field ob- The standard deviations of the averages are given 
by the blue solid circles or blue solid diamonds, respectively. 
Here we take Bq = 0. 



suppressed. 

InRef.0, we have shown that a constant external mag- 
netic field makes the second-order effective interaction 
mediated by an even-size bus anisotropic. Here we find 
that local random variations in the magnetic field, 5 Bi, 
also induce anisotropy, even if Bq = and (8Bi) cn = 0, as 
shown in panel (c) of Fig. 1121 This anisotropy is weaker 
compared to the case of an odd bus in a finite Bq as 
shown in Fig. [9l because the average of the local field 
here vanishes (Bq = 0), so that the spatial isotropy is 
not broken as completely as in the case of Fig. [5] Note 
that we would normally expect an even-size bus to be 
operated at zero external magnetic field Bq = 0, unless 
the anisotropy of the effective coupling is desired. 



try also makes the ground Zeeman energy gap of the odd- 
size bus robust against small fluctuations in the exchange 
couplings. However, randomness in the exchange cou- 
plings does cause the eigenenergies and the eigenstates 
to vary, which in turn leads to randomness in the mag- 
nitudes of the effective couplings (both qubit-bus and 
qubit-qubit). 

An external magnetic field, whether uniform or ran- 
dom, does break the isotropy of the Heisenberg spin bus, 
and leads to anisotropy in the qubit-bus and qubit-qubit 
effective couplings. A locally random magnetic field also 
lifts the ground state degeneracy of an odd-size bus, even 
when the average applied field vanishes. The local ran- 
domness also gives rise to the effect that the total mag- 
netic moment of the odd size bus in the ground state 
may be antiparallel to the direction of the applied mag- 
netic fields, when a single spin in the ground state would 
have been parallel to the magnetic field. Even-size buses 
are somewhat more robust against local random mag- 
netic fields, since their ground state is non-magnetic. 

We have performed ensemble calculations for the cou- 
pled qubit-bus systems we have considered, where the 
standard deviation of an ensemble averaged quantity 
(such as the qubit-bus and qubit-qubit effective cou- 
plings) represents the sensitivity of this quantity to the 
particular parameter randomness. Thus our results have 
clear implications not only for situations where static pa- 
rameter randomness is present, but also for dynamical 
noise in the exchange coupling or the external field. 

To give context to this paper, in the Appendices we 
provide a comprehensive overview of even and odd spin 
chains as quantum data buses. In particular, we derive 
the first- and second-order effective couplings using the 
projection operator method. We explore the low-energy 
spectra of buses coupled to zero, one or two qubits, 
from which we derive the first- and second-order effec- 
tive Hamiltonians of the qubit-bus system. We also prove 
that random exchange couplings do not lift the ground 
state degeneracy of an odd-size bus. Finally, we present 
a study of the scaling properties of the bus, for up to 20 
nodes. 
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IV. CONCLUSIONS 

We have performed a comprehensive study of the 
effects of local randomness in the exchange couplings 
and the external magnetic field on the capabilities of a 
strongly coupled Heisenberg spin bus. 

We find that the random exchange couplings preserve 
the isotropic symmetry (in the qubit-bus and qubit-qubit 
couplings) of the bare Heisenberg coupling. This symme- 



Appendix A: Derivation of the spin bus effective 
Hamiltonians by projection method 

In this Appendix we derive the effective low-energy 
Hamiltonian of the full Hamiltonian (|2.ip . The weak 
qubit-bus couplings J a ,i/ Jq <C 1 are used as perturbation 
parameters. The total Hamiltonian H = H + Hi can be 
rewritten as an unperturbed Hamiltonian H = Hq+Hq 
and a perturbation H\ = Hqc- Since we are interested 
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in the low-energy limit, we define a projection operator 
P onto the subspace Ho spanned by the tensor products 
|$fc) of the ground state(s) of the free bus Hamiltonian 
He and the eigenstates of the free qubit Hamiltonian Hq 

p= i**x**i- (Ai) 

k£Ho 

Here the Zeeman energy of Hq is assumed to be small 
compared to the gap of He, so that the energy of the 
subspace Ho is equal or nearly equal to the ground energy 
Eq of He- In this limit, we can set Hq = 0. The finite- 
field induced anisotropy is addressed elsewhere £22 The 
effective Hamiltonian acting on the subspace Ho is then 
given b y 16 ' 17 

H cB = PHP + PHQ- L— QHP, (A2) 

tiQ — LjtL Q 

where Q = I — P projects onto the sub-Hilbert-space 
orthogonal to Ho- The effective Hamiltonian to second 
order in Hi is given by 

H$ = PHP + PH X Q - 1 - QHyP . (A3) 
hiQ — Ho 

The derivation of the explicit form of the effective Hamil- 
tonian (|A3|) requires detailed information on P, which 
consists of the structure and spectrum of the ground 
manifold of He- 

The Heisenberg spin chain Hamiltonian He, given by 
Eq. (|2.2[) , is exactly albeit only partially solvable with 
the Bethe ansatz.— The z component S z of the total 
spin S = Si commutes with the Hamiltonian (|2.2[) . 
[H ,S Z ] = 0, so that the energy eigenstates can be la- 
beled by \n\ m z ) with the energy level n and the magnetic 
quantum number m z , i.e., the eigenvalues of S z . How- 
ever, the general analytic expressions of the eigenstates 
|n; m z ) are not available. Traditionally, for bulk systems, 
the periodic boundary condition and an even number N 
(N — > oo in thermodynamic limit) are assumed. For fi- 
nite size chains, however, the eigenstates are dependent 
on both the boundary condition and the even- odd parity 
of size N, as shown in Figs. [13] and [15] Thus, the ef- 
fective Hamiltonians for the bus-qubit system are differ- 
ent depending on the parity of the bus£ Here we give 
a more detailed description of the derivation of the ef- 
fective Hamiltonians. Note that Ref. [lj^ demonstrated 
experimentally the even-odd parity effect of a spin chain, 
by assembling chains of 1 to 10 Mn atoms on a metal- 
lic surface and measuring the parity dependent tunneling 
currents. 



1. Effective Hamiltonians with an Odd-Size Bus 

An odd-size antiferromagnetic chain has an odd num- 
ber of spins, so that the ground state should have one un- 
compensated spin. For example, for N = 3, the classical 
antiferromagnetic spin configurations are "up-down-up" 



or "down-up-down" . The exact degenerate quantum me- 
chanical ground states of the odd-size chain with N = 3 
around Bo — are given by 

|0; +i) c = -L (|001) - 2|010) + |100» , (A4a) 
|0; -i) c = -L (|H0) - 2|101) + |011)) , (A4b) 

where |0) and |1) on the right hand side represent the spin 
up and down states of a single spin, respectively. One can 
see that the basis states corresponding to classical con- 
figurations, |010) or 1 101), are most probable, although 
quantum corrections are already sizable. For longer 
chains, the amplitude of the classical antiferromagnetic 
configuration continues to decrease, while quantum fluc- 
tuation contributions increase. Although a three-node 
chain is small, the analytic solution can serve as a start- 
ing point for understanding how an odd-size chain acts 
as an effective single spin. 

To study longer spin chains, we numerically solve the 
eigenvalues and eigenvectors of the Hamiltonian (|2.2|) 
with LAPACK.— Figure [T3l (a) plots a few lowest energy 
levels of a spin chain with N = 7 nodes as a function 
of the external magnetic field B - When B = 0, the 
odd-size chain has two doubly degenerate ground states, 
with total magnetic quantum number S z = ±1/2, and 
denoted by |0; ^)c and |0; — \)c- An external magnetic 
field Bq splits the two degenerate ground states by the 
Zeeman energy (for small fields), labeled by Aqi. Thus, 
the odd-size spin chain can be considered as an effective 
single spin at small Bo field and in the low energy limit, 
as shown in Fig. Q2] (b). 

Two important energy scales for an odd-size chain are 
the energy gaps A i and A 12 . The energy gap A 12 gives 
the splitting between the ground and excited manifolds. 
The fact that Aoi is given by the Zeeman splitting at 
low fields signifies that the doublet can be considered an 
effective spin-1/2. Figure [14] shows how the description 
of the odd-size chain as an effective single spin is lim- 
ited by its size N. While the Zeeman splitting Am is 
independent of the size N of the spin chain around zero 
magnetic field, the range of the magnetic field -Bo, within 
which the effective spin picture is valid, does depend on 
N . The crossover behavior in Aoi occurs at smaller Bo 
for larger N because the manifold splitting, A12 ~ 1/iV, 
decreases with N. 

Now that we have established the fact that an odd- 
size spin chain can be considered as an effective spin-1/2 
at low energy, we can construct an effective qubit-bus 
Hamiltonian within the Hilbert space spanned by the bus 
ground doublet states and the qubit states. With the 
doubly degenerate ground states |0;±i) of the odd-size 
bus, the projection operator P is written as 

P = Ia,b ® ( |0; i)(0; ±| + |0; -±)(0; -i| ) c (A5) 

The effective Hamiltonian to the first order in Hqc is 

H$ = PHP = PH P + PH X P , (A6) 
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FIG. 13: (Color online) (a) A few lowest energy levels of an 
odd-size chain with TV = 7 as a function of Bq. The two 
lowest levels around So form up and down states of an effec- 
tive single spin. Two gaps, Aoi and A12 are indicated, (b) 
A magnified plot of the two lowest levels, offset against the 
ground state energy Eo/Jo — —2.83624 at Bo = 0, shows a 
Zeeman splitting of a single spin 1/2. 



where the effective coupling between qubit a and the 
spin-bus C is given by 



T^ 1 ) — J 777 

rrii = (0; ±|er iz |0; \) = 
= (0;i|a«|0;-i; 



-(0;-i|a«|0;-i) 



(A8a) 
(A8b) 



where J a A is the bare coupling between the ith spin of 
the chain and the external qubit a, and mj is the dimen- 
sionless local magnetic moment at the site i of the chain 
in the ground state. In the case of N = 3, Equation (|A4|) 
gives the local magnetic moments, mi = 2/3, 777,2 = 
— 1/3, 7773 = 2/3, with mi + 777 2 + 7Ti 3 = 1. The effective 
Hamiltonian (|A7[) shows again that an odd-size chain 
acts as an effective spin- 1/2 particle that couples to the 
external qubits A and B, as illustrated in Fig. [I] 

The- effective Hamiltonian to second-order in Hi. which 
is needed for longer-time operations, is given by£ 



hi 



(A9) 



where the RKKY-like second-order coupling J^ 2 ) igjl&lL 



J. 



(2) 



JaaJb, 



\a itl \n)(n\(T m \Q) 

Eq — E n 



(A10) 
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FIG. 14: (Color online) Two energy gaps, Aoi and A12 
measured in units of Jo are plotted as a function of Bo for 



where PHqP gives rise to a constant energy shift. In an 
external magnetic field, it also gives the Zeeman interac- 
tion of the effective spin- 1/2 system of the spin bus or 
qubits, g[iB B • Sc or g^s B • S Q with the external field. 
With Bq — 0, the first-order effective Hamiltonian^ is 
only given by PH\P, and takes the form 



H 



(i) 

eff 



J A A a A 



S c + S B 



(A7) 



Here E n are the energy levels and \ri) are the eigen- 
states of He omitting the magnetic quantum number. 
The summation with a prime indicates the exclusion of 
the ground states. While the exact calculation of J- J re- 
quires complete information of the eigenvalues and eigen- 
states of the chain, an approximate form can be obtained 
with only the ground state and the excitation gap. Using 
the closure relation X) l n )( n l = 1) where n sums through 
all the isolated-chain eigenstates, we obtain 



J$/(JA,iJB,j) 

1 -((0\o- iz 



2Ai 



a JZ \0)-(0\a lz a jz \0) 



(All) 



where |0) refers to |0;±i)p. In other words, the second- 

order coupling J>J can be approximated by the difference 
between the spin-spin correlation function and a product 
of local magnetic moments of the ground state of the 

odd-size chain. The competition between these two terms 

(2) 

leads to decaying oscillations in J>J. Note that a phase 
slip occurs when the qubit separation reaches a certain 
range, as discussed in Ref. H 



2. The Effective Hamiltonian with an Even-Size 
Bus 

In an even-size chain with antiferromagnetic couplings, 
the spins are completely compensated, so that the chain 
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has a non-degenerate ground state with both zero total 
magnetic moment, \n = 0; S z = 0)c, and zero local mag- 
netic moment, (0; 0|o" p i |0; 0) = 0. These properties of 
the even-size chain can be illustrated with the simplest 
case of N = 2, when the ground state is the singlet state 
\n = 0;S z = 0) = ^(|0, 1) - |1, 0)). In Fig. M we plot 
the lowest energy levels of an even-size chain with N = 8 
as a function of the external magnetic field Bq. For small 
Bo, the non-degenerate ground state with S z = is sep- 
arated from the excited states by the gap Aoi. Similar 
to the odd-size chain, this gap and the field range before 
level crossing is limited by the size of the chain, as shown 
in Fig. HH 
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FIG. 15: (Color online) The four lowest energy levels of an 
even-size chain with N = 8 as a function of Bo. The two 
lowest energy gaps, Aoi and A12 in units of Jo are indicated. 
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FIG. 16: (Color online) The two lowest energy gaps, Aoi and 
A12, of even-size chains with N = 4, 6, 8, 10, 12, as functions 
of B . 

Protected by the Aoi gap, the low-energy Hilbert space 
we are interested in is spanned by the ground state of the 
even-size chain and the qubit basis states. The projection 
operator takes the form P = l AB <g> (|0; 0)(0; 0|) c . The 
first order term in Eq. (JX3) is PHP = PH Q P + PHiP = 
EqP, which is just a constant energy shift. In other 
words, the first order perturbation term for an even- 
size chain induces no qubit-bus couplings, since the bus 
ground state by itself cannot have any dynamics. The 



second-order perturbation term is obtained in the same 
way as in the case of an odd-size chain. The effective 
Hamiltonian to second order in H\ is 6 



H 



(2) 
elf 



J 



(A12) 



(2) 

where the RKKY-like coupling J\j has the same form as 
Eq. (jAlOj) . although the prime here would exclude only 
a single ground state. Similar to the case of the odd-size 
chain, J>j can be approximated as 



1 



2A 



01 



(0;0\a iz a jz \0;0), (A13) 



where Aoi is the energy gap between the ground state 
and the first excited state, and (0; 0|<7i Z crj Z |0; 0) is the 
spin-spin correlation function of the ground state of the 
chain. There is no contribution from local magnetic mo- 
ments here since they vanish in the ground state of an 
even-size bus. 



Appendix B: Absence of effect on ground state 
splitting of an odd-size bus by random exchange 
couplings 

Here we prove that for an odd-size spin chain, the en- 
ergy splitting A i of the ground doublet states at low 
uniform magnetic fields is independent of the randomness 
in the inter- node exchange coupling J^. With E>i = Bq, 
the Hamiltonian (|3.ip is 



H C = 4 0) 



V. 



(Bl) 



where Hq is the unperturbed Hamiltonian (|2 . 2[> . and 
the perturbation V is given by 



V 



N-l 

E 

i=l 



6Ji Si ■ Sj+i 



(B2) 



The two lowest eigenstates of the unperturbed Hamil- 
tonian H { ° } are |0; -f) = |0^) and |1; ±) = |l( )), when 
the magnetic field is applied in the positive z direction, 
and the corresponding eigenenergies are and ^ 0) . 
For the perturbed Hamiltonian He, we denote the two 
lowest eigenstates by |0) and |1), respectively, and the 
corresponding eigenvalues by Eq and E\. By definition, 
the energy gap Aqi of He can be expressed as 



Aoi = E-l - E = + (5E 1 - 5Eq) , 



(B3) 



x ni 



E^ ' — Eq' 1 is the lowest energy gap of the 



h(0) 



unperturbed Hamiltonian \ which is a Zeeman gap. 
The energy shifts of the two lowest levels, 5Eq = Eq—E^ 
and 5 Ei = E± — e[°\ caused by the perturbation V, i.e., 
fluctuations SJi, are given by^i 

(0<°>|V|0) 



se 



(o(°)|o) 



SEi = 



(l (0) l^|l) 
(l(o)|l) 



(B4) 
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Their difference 8E\ — SEq is thus 

(1 (0) |F|1) (o ( °)|y|o) 



5E 1 - SE = 

N-l 



!=1 



(1(0)|1) 

l (0) |s»-s J+1 |l) 



awn) 



(0(°) |0) 

(0W|s. t -s m |0) 
(0(°)|0) 



(B5a) 



(B5b) 



Since [Hc,S z ] = 0, the two lowest states, |0) and |1), 
(also |0^ ') and ll^ -*)) of an odd-size bus are spin- flipped 
states of each other, so that the bracket part in the above 
equation vanishes, which leads to SEi — SEq = 0. This 
means that the two lowest energy states fluctuate to- 
gether, as shown in Fig. [5] (a), so that their difference, 
i.e., the Zeeman energy gap, does not change. 

We have thus proven that the Zeeman energy splitting 
A i of the ground doublet states is invariant over random 
exchange couplings Jj. Consequently, the odd-size chain 
with random exchange couplings can still be regarded as 
an effective single spin 1/2 in the low energy limit. 



Appendix C: Scaling Properties of Spin Buses 

In this Appendix we discuss the scaling properties of a 
spin chain as a quantum data bus. While the Heisenberg 
spin-1/2 chain is exactly solvable with Bethe ansatz^ 
only partial information about the ground state and the 
elementary excitations are available. Various numerical 
approaches have been applied to this system since Bon- 
ner and Fisher's pioneering work^ Although there has 
been tremendous advances in computational power, the 
exact diagonalization method can only handle a spin-1/2 
system with sizes of up to N ~ 40, depending on the 
number of eigenvalues and eigenstates to be calculated. 
Indeed, one could consider such limitations as one of the 
motivations for building a quantum computer. Below we 
present our results for spin chains with N up to 20. 

Calculation of the first-order qubit-bus effective cou- 
pling of Eq. (|A8p needs only knowledge on the ground 
state of an odd-size chain. On the other hand, calcu- 
lation of the second-order coupling, Eq. (|A10[) . requires 
the full knowledge of the eigenvalues and eigenstates of 
the spin chain. We have solved the full spectrum of spin 
chains with sizes of up to N = 14 on a personal com- 
puter using LAPACK— , and obtained a few lowest en- 
ergy eigenvalues and eigenstates for spin chains with sizes 
up to N = 20 using PRIMMED 

Figure [17] shows our results on the ground state en- 
ergy Eq per bond, the gap between the ground and the 
excited states, and the end-point local magnetic moment 
of the odd-size chain for finite chains of size N = 2 to 
20. As shown in Fig. [IT] (a), the ground state energy 



per bond approaches the analytic value of Eo/(NbJo) = 
i— hi 2 ps —0.4431 in the thermodynamics limit, obtained 
from the Bethe ansatz. Here the number 2V& of bonds is 
given by JVj = N — 1 for open chains and Nf, = N for 
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FIG. 17: (Color online) (a) The ground state energy Eo /Jo 
per bond of the open chains (open red circle) and rings (solid 
blue circle), (b) the gap A between the ground state and the 
excited state, (c) the local magnetic moment tun at the end 
of the open-chain as a function of the size N, and (d) a log- 
log plot of mjv for odd-size chains. The dashed line in panel 
(a) indicates E /(J N) = | — In 2 ps —0.4431 obtained from 
the Bethe ansata^. The dashed line in panel (b) represents 

The gray lines in panels (c) and 
N fit. 



rings. The ground state energy oscillates, depending on 
the even-odd parity, as the size of the chain N increases. 
However, this finite-size effect diminishes in the large- N 
limit. Figure [17] (b) plots the ground state energy gap A 
as a function of the size of the chain N, with A = Aoi for 
even-size chains and A = A12 for odd-size chains. The 
numerical result follows the well-known analytical esti- 
mate A ~ ir 2 Jo/2N as the size N increases^ Although 
the ground energy E , the spin-spin correlation function 
for the ground state (sj • Sj), and the ground energy gap 
A are well known ; 26 ' 27 the scaling property of the lo- 
cal magnetic moment of the odd-size chain is less well 
understood.— Figure [ijj (c) plots the dependence of the 
end-site local magnetic moment on the size N. Our nu- 
merical data points to a rn/v oc 1/y/N dependence (pos- 
sibly slightly faster), as indicated in panel (d) of Fig. [17] 
Further work is still needed to clarify this ./V-dependence. 
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